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ABSTRACT 

Context. The study of pulsation in Pre-Main-Sequence intermediate-mass stars represents an important tool for de- 
riving information on fundamental stellar parameters and internal structure, as well as for testing current theoretical 
models. Interest in this class of variable stars has significantly increased during the last decade and about 30 members 
are presently known in the literature. 

Aims. We have constructed the frequency spectrum of the oscillations in V346 Ori. We apply asteroseismic tools to 
these data to estimate the intrinsic parameters (mass, luminosity, effective temperature) of V346 Ori and to obtain 
information on its internal structure. 

Methods. CCD time series photometry in the Johnson V filter has been obtained for a total of 145.7 h of observations 
distributed over 36 nights. The resulting light curves have been subjected to a detailed frequency analysis using updated 
numerical techniques. Photometric and spectroscopic data have also been acquired to determine reliable estimates of 
the stellar properties. 

Results. We have identified 13 oscillation frequencies, 6 of which with higher significance. These have been compared 
with the predictions of non-radial adiabatic models. The resulting best fit model has a mass of 2.1±0.2 Mq, luminosity 
logL/I/Q = 1.371q5^3, ^'id effective temperature 7300±200 K. These values are marginally consistent with the associa- 
tion of V346 Ori to Orion OB la. Alternatively, V346 Ori could be placed at a slightly larger distance than previously 
estimated. 

Key words, stars: variables: S Scuti - stars: oscillations - stars: pre-main sequence - stars: fundamental parameters - 
stars: individual V346 Ori 



1. Introduction 

Pre-Main-Sequence (PMS) stars with masses larger than 
~1.5 M0 that show variable emis sion lines are known as 
Herbig Ae/Be stars ()Herbiel Il960f ). Observationally, they 
are found within star forming regions and show variable 
emission lines (especially Ha) and strong infrared excess 
caused by the presence of circumstellar material. In addi- 
tion, Herbig Ae/Be stars are characterized by photometric 
and spectroscopic variability on time scales of minutes to 
years, mainly due to photospheric activity and interaction 



Send offprint requests to: V. Ripepi 

* Based on data collected with the REM 0.6-m, Loiano 1.5-m, 
SARA 0.9-m and SNO 1.5-m telescopes 



with the circumstellar environment (e.g.. lGahm et al]|1995l : 
iBohm fc Catalalll995l) . 

Although recent theoretical work has i mprove d our un 



derst a nding of PMS evolution (see e.g. iPalla fc Stable: 



11993t iDAntona fc Mazzitellil 11994 ISwenson et all 1199 
there remain differences in the models owing to uncertainty 
in the treatment of convection, of the input physics and 
the zero-point of the stellar ages. It is therefore desirable 
to find independent ways to constrain PMS evolutionary 
tracks and, in turn, the internal structure of intermediate- 
mass stars which are subject to dramatic changes during 
the short-lived PMS contraction phase. 

Asteroseismology of Herbig Ae/Be stars can in principle 
test PMS models by probing their interiors. It is now well 
established that these stars cross the pulsation instability 
strip of more evolved stars during contraction toward the 
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main sequence, co mprising a class of variable stars called 
PMS S Scuti stars ("Kurtz fc Marang'l995 HMarconi fc Pallal 
[1998 : Catala 2003; Ripcpi ct al. 2006a). 



The first theoretical investigation of the PMS instability 
strip based on nonlinear convective hydrodynamical mod- 
els was carried out by Marconi & Palla (1998, hereinafter 
MP98) who calculated its topology for the fi rst three ra- 
dial m odes. A subsequent theoretical work bv lSuran et al.l 
(|200lh made a comparative study of the seismology of a 
1.8 Mq PMS and post-MS star. These authors found that 
the unstable frequency range is roughly the same for PMS 
and post-MS stars, but that some non-radial (g) modes are 
very sensitive to the dee p internal structure. More recently, 
iGrigahcene et all (j2006l ) have produced a theoretical insta- 
bility strip for PMS stars for the first seven radial modes. 
Since the work by MP98, many new PMS S Scuti candi- 
dates have been observed and the current census of known 
or s uspected candidates includes about 30 stars (see e.g. 
iRipep i et al. 2006a; Zwintz 2008). Ho wever, only a few of 
these have been s t udied in detail (e.g. iRipepi et al.ll2003t 
iBohm et al.l 12004 IRipepi et aLll2006lj ). Thus, the overaU 
properties of this class of variable stars are still poorly de- 
termined. 

In order to perform reliable asteroseismological analy- 
sis, it is critical to obtain data more accurate than that 
available in the literature. In particular, it is necessary to 
measure as many excited oscillations as possible with high 
precision on frequency and, ideally, absence of aliases. Prior 
to space data which will become progressively available 
from the MOST, COROT and possibly KEPLER satel- 
lites, the only way we can satisfy these requirements is 
to carry out multisite observations on promising astero- 
seismological targets, i.e. stars which are multiperiodic. 
Starting in 2003, our group has applie d the mu l tisite tech- 
nique t o two targets, V35 1 Ori ( Rip epi et al.l I2003D and 
IP Per (|Ripepi et al.ll2006bf ). The observational results ob- 
tained on these stars have been recently int erpreted using 
non-r adial adiabatic pulsation theory (see iRuoppo et al.l 
|2007() . As a continuation of our ongoing project, here we 
concentrate on the Herbig Ae star V346 Ori. This object 
is cla ssified as a member of the Orion OBla subassocia- 
tion (jHernandez et al.ll200 6l). or of its subgroup centered 



on 25 Ori (jBriceno et al.l l2007). It w as originally su s pecte d 
to be variable on short time s c ale by [Marconi et all (|2000D . 
Subsequently, iPinheiro et al.l (|2003[ ). on the basis of only 
two observing nights, confirmed the S Scuti nature of the 
light variation by measuring two highly significant frequen- 
cies of pulsation at /i=35.2±1.3 c/d and /2=22.6±1.5 c/d 
and two less significant frequencies at /3=45.7±1.3 c/d and 
/4=18.7±1.3 c/d. The measured amplitudes of the oscilla- 
tions in the V band were Ai =2.9, A2 —2.8, A3 —1.2, and 
A4 =1.4 mmag. To improve on these initial results, we car- 
ried out: 1) an extensive photometric multisite campaign 
to enlarge the number of observed mode of pulsation; and 
2) a spectroscopic analysis to determine the values of log g, 
Tgs and v sin i of V346 Ori; 3) a comparison with nonlinear 
pulsation models to constrain the position of V346 Ori in 
the HR diagram. 

The structure of the paper is as follows: in Sect. 2 we 
describe the observations and data reduction; in Sect. 3 we 
deal with the frequency analysis; in Sect. 4 we report the 
results of the spectroscopic analysis; in Sect. 5 we compare 
the measured frequencies with model predictions and derive 
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Fig. 1. Fourier analysis of the differential photometry 
(TYC-101-317-1)-(HD 287842). See the text 

an estimate of the parameters of V346 Ori; and in Sect. 6 
we give some final remarks. 

2. Observations and data reduction 

2.1. Time series photometry 

Time series observations were obtained using four differ- 
ent telescopes equipped with CCD devices during 36 nights 
between 2005, October and 2006, January. A complete de- 
scription of the observations is given in Table lA.ll (only 
available electronically). As comparison star, we observed 
HD 287842 (J2000: 05'' 24" 24.2P -f01° 4l' 33.3"; V=10.3, 
A7) . This star was chosen because it is only ~5 away from 
V346 Ori and within the field of view of all the instruments 
used. Also, its magnitude and color are not too different 
from those of V346 Ori. The constancy of HD 287842 has 
been checked with respect to other (fainter) stars present 
in the CCD frames and in particular with TYC-101-317- 
1, a bright star (V~11.5 mag) which is present in the 
Sierra Nevada data (five nights, see Table lA.ll and below 
for the data analysis) . The Fourier analysis of the differen- 
tial photometrjO (TYC-101-317-1)-(HD 287842) is plotted 
in Fig. [TJ No significant frequency peak is present in the 
periodogram. The observations were carried out only with 
the V filter to have a time sampling sufficiently fast to sam- 
ple properly the expected rapid pulsations. In the following 
we describe the main features of the instrumentation used 
for this work and give some details concerning the data 
reduction and photometric procedures. 

— Loiancl3: the observations were carried out using the 
1.5-m telescope equipped with the BFOSC instrument. 
The CCD was a EEV 1340x1300 pixels with individual 
size of 0.58 , for a total field of view of 13 x 13 . The 
data have been reduced following the usual procedures 
(de-biasing, flat-fielding) and using standard IRAF rou- 
tines. The aperture photometry has been carried out 
using routines written in the MIDAS environment. 

— REIVQ: the data were collected using the 0.60-m au- 
tomated telescope at La Silla (Chile). Images were ob- 
tained with the ROSS spectrograph equipped with a 
commercial Apogee AP47 camera hosting a Marconi 47- 
10 IKxlK 13 nm pitch CCD, which covers 9.54' x 9.54' 
with a scale of 0.56 /pixel. For the data reduction and 



^ Note that small nightly variations were removed by sub- 
tracting for each night the average photometry. 
^ http:/ /www. bo. astro.it/loiano/index. htm 
^ http://www.rem.inaf.it 
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Fig. 2. Light curves for V346 Ori of the first 20 nights. Note that SV = Vvar ~ Vcomp, where VAR refers to V346 Ori 
and COMP is the comparison star HD 287842. In aU panels, the solid line displays the fit to the data with all the 
significant frequencies listed in Table ETT] fsee text). 



photometry we followed the same recipes as for the 
Loiano observations. 
— SARA|3 (Southeastern Association for Research in 
Astronomy): we used the 0.9-m automated telescope 
at Kitt Peak (USA). The images were collected in the 
V filter using an Apogee AP7p camera with a back- 
illuminated SITe SIA 502AB 512x512 pixel CCD. The 



pixels are 24 microns square, corresponding to 0.73 at 
the telescope focal plane scale, giving a field of view of 
6.2 X 6.2 . Sky fiats, dark and bias exposures were taken 
every night. The data were calibrated and reduced using 
standard IRAF routines. 
— SNO0 (Sierra Nevada Observatory): the data was ob- 
tained using the 1.5-m telescope equipped with a 2k x 2k 



* http://saraobservatory.org 



^ http://www.osn.iaa.es 
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Fig. 3. As in Fig. 1, but for the remaining 16 nights. 



CCD, with a pixel size of 0.23 for a total field of view 
of 7.92 X 7.92 . Observations were carried out in the 
V filter with exposures of about 7 s (depending on the 
quality of the night), using a 2x2 binning. The data were 
reduced in a standard way by correcting the frames from 
bias and flat-field. Four different apertures were tested 
with radii of 4, 8, 12 and 16 pixels. The 12 pixel aperture 
(corresponding to about 5 ) was chosen as the optimal 
for our observations. The sky inner radius was fixed to 
23 pixels and the annulus width to 3 pixels. 

In total, we collected 145.7 hr of observations during 36 
nights, distributed as follows: 18.7 hr with Loiano; 94.3 hr 
with REM ; 11.2 hr with SARA; 21.5 hr with SNO. 

To prepare the data for Fourier analysis, we have first 
detrended the data to a common average zero value. This is 
needed for the following reasons: 1) to avoid problems due 
to differences in the instruments-|-filters used; 2) to pre- 
vent problems with zero point differences between different 
datasets; 3) to eliminate the long term photometric varia- 
tions due to the presence of circumstellar material around 



the star, a well known occurrence for Herbig Ae stars such 
as V346 Ori. The detrending procedure consisted in remov- 
ing a constant from the photometry for the majority of the 
nights. In the remaining cases, for nights with a duration 
of the observations larger than approximately 5 hours, we 
used a linear or quadratic regression to dctrend the data. As 
a consequence of this procedure, we are not able to investi- 
gate frequencies lower than about 4.5-5 c/d. The resulting 
light curve of V346 Ori is shown in Fig. [2] and [2] Here, the 
solid lines represent the least square fit to the data based 
on the frequency analysis discussed in Section 3. 

2.2. Stromgren-Crawford photometry 

A few uvhyP points were also collected for V346 Ori 
in the Stromgren-Crawford photometric system for cal- 
ibrating purposes using the six-channel uvby/3 spectro- 
graph attached to the 0.9-m telescope at Sierra Nevada 
Observatory. During these measurements, three bright 
stars C1=HD 34888, C2=HD34745 and C3=HD 36525 
were used as check stars, and instrumental magnitude 
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Table 1. uvhy(3 indices obtained for V346 Ori and compar- 
ison stars. The pairs below the star names are the number 
of points collected for each object in uv hy and (3, respec- 
tively. The values given bv loIsCTl (fl996h are listed in the 
bottom part. 
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Fig. 4. Spectral window (in power) 



differences were obtained relative to Cl. To transform 
these instrumental differences into the standard uvby(3 sys- 
tem a set of 19 standards stars wer e selected from the 
list of ICrawford fc Mandal (|1966D and lCrawford fc Barnes! 
(fT970l) . 

Next, the absolute standard uvhyfi indices for V346 Ori 
and check stars were obta i ned f ollowing the method de- 
scribed in [Rodriguez et all (j2003[ ). using Cl to C3 as zero- 
points. The results are listed in Table [1] together with the 
values reported in the literature for the check stars. The 
error bars in this table mean standard deviations of mag- 
nitude differences relative to Cl. As seen, our results are in 
very good agreement with t he va lues found in the homo- 
geneous catalogue of Olscn (1996). Simila r results can be 
find in the lists of .Olsen (1983^ and .Hauck fc MermilliodI 

diaii). 

These indices, together with suitable calibrations avail- 
able in the literature for uvbyf3, can be used to estimate the 
physical parameters of V346 Ori using the procedure de- 
scribed in [Rodriguez et all (|2001h . In particular, a value of 
(5mi=0.012 is obtained for the metallicity parameter which 
leads to nearly so l ar ab undances of [M/H]=— 0.05(±0.1) 
using the ISmallevI (|l993i' ) relati on for metal abun d ances . 
Moreover, using the grids of ISmallev fc Kupkal (|1997[ ) 
for solar abundances, values of Teff=7340(±150) K and 
log5=4.0 (±0.1) are also found which are in very good 
agreement with those obtained from spectroscopy and dis- 
cussed in Section 4. 

2.3. Spectroscopic data 

An echelle spectrum for V346 Ori was obtained with the 
BFOSC instrument at the Loiano Telescope during the 
night of January 31, 2007. We used the grism #9 with 
the cross disperser #10 in order to improve the efhciency 
towards the blue. The resulting spectrum includes 13^ or- 
ders and covers the wavelength interval from ^3800 A to 
~10000 A. Data reduction has been performed using the 
IRAF package ECHELLE. The resolution of our data was 
R « 2500. The spectrum of V346 Ori is shown and discussed 
in Section 4. 



3. Frequency analysis 

The pulsation frequenc y analysis has been ca. rried out with 
the periodOA package ()Lenz fc Bregerl l2005t ) that adopts 
both Fourier and least-squares algorithms and permits si- 
multaneously fitting of multiple sinusoidal variations, thus 
it does not rely on sequential prewhitening. In addition, 
periodOA offers the possibility to weight the data during 
the least-squares procedure (see below). 

Before analysing the time series, we first note that 
the photometric quality of the data varies among different 
nights and sometimes during the same night. To account 
for this and to include the measured photometric errors 
based on Poisson statistics which often are underestimated 
with respect to the real scatter of the data, we have used 
point-to-point weights for the frequency analysis. Thus, we 
adopted the weighting scheme in periodOA called "deviation 
weight" . Defining residuals as the absolute value of the 
difference between observed and calculated phase points, 
this procedure assigns the weights as follows: weight=l 
for residual < cutoff; weigh.t=[{cutof f /residual)^] for 
residual > cutoff. The choice of the cutoff value is crit- 
ical to obtaining reliable results. A general rule is to avoid 
using a cutoff to o small compa red to the amplitude of the 
oscillations (e.g. iHandlerj |2003| ). We found that a cutoff 
of 7.5 mmag complies with these requirements and verified 
that lower values give unrealistic results, whereas higher 
values (e.g cutof f=\ mmag) do not produce results sig- 
nificantly different than the 7.5 mmag case. As described 
below, we performed the frequency analysis in two ways: 
with and without weights. 

The frequency analysis requires first the discussion of 
the Spectral Window (SW) which is shown in Fig. [4l 
Unfortunately, due to the lack of an observing site in Asia, 
we were not able to completely cancel out the 1 c/d alias. 
Nevertheless, the SW of our data appears of good quality in 
both the weight and no-weight case, as the 1 c/d is reduced 
to about 50% (in power), and the 2 c/d alias is almost 
negligible. After this, the procedure consists of the follow- 
ing steps: i) calculation of Fourier transform and identifi- 
cation of the highest amplitude peak; ii) least-square fit to 
the data with this frequency, improving simultaneously fre- 
quency, amplitude and phase; iii) check on the significance 



6 



S. Bernabei et al.: Multisite photometry of the pulsating Herbig Ae star V346 Ori 



Table 2. Frequencies and amplitudes for V346 Ori calculated with periodOA for the two cases, with and without weights 
(see text). The estimated error on frequency is ~ ±0.02 c/d (0.23 fiHz). For the no-weight case the frequencies are 
listed in the same order of the weighted mode and not by decreasing amplitude. The horizontal line separates the high 
significant frequencies from the less significant ones. 
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of the extracted frequency; iv) calculation of Fourier trans- 
form of the residuals and identification of the new highest 
amplitude peak; v) least-square fit to the data using simul- 
taneously the derived frequencies; vi) iteration of the previ- 
ous steps for the identification of all significant frequencies. 
To determine whether to accept or reject a frequency, we 
adopted the widely used requiren ient of a signal to noise 
(S/N) amplitude ratio equal to 4 (jBreger e t al."l993 ), cor- 
respo nding to a value of about 12.6 in power (Brcge r et al.l 
I2005D . In order to calculate the noise level, we used again 
periodOA to obtain an average value for the noise itself in 
bins of 10 c/d width, moving along the spectrum. We ver- 
ified that changing the size of the bins up to 50% does 
not alter the S/N level significantly, especially for the low- 
amplitude frequencies close to the significance limit. 

The result of the frequency analysis is shown in Fig. [SI 
while the derived frequencies are listed in Table [2] for the 
weight and no-weight case. The difference in the two cases 
is significant, since the number of derived frequencies varies 
from 13 in the weighted case to only 6 in the non- weighted 
case. However, note that five of the six derived frequencies 
coincide in both cases, while /g in the non-weighted model 
is clearly the -1 c/d alias of /g in the weighted mode. As 
seen in Table [21 the use of weights has the effect of reduc- 
ing the amplitude peaks. More significantly, the noise is 
also reduced, allowing more frequencies to exceed the Aa 
threshold albeit only marginally (see fr to /13 in Table [2]). 
Given the poor significance of the frequencies derived with 
our weighting scheme, in the comparison with the models 
we shall use only the first six frequencies in Table [31 The 
remaining ones are in any case useful to verify if the best-fit 
model is also able to reproduce them. 

Finally, we have checked if some frequency could simply 
be a linear combination of others and verified that all the 
detected periods are independent of each other. 

Regarding the errors on the extracted frequencies, we 
adopted the classical method of calculating the FWHM of 
the main lobe in t he SW, obtaining af ~0.02 c/d. (see 
[Alvarez et al.|[T998l for a discussion on the subject). 



We ca n now compare the present results with those ob- 
tained bv lPinheiro et al.l (|2003). As a result, the two signif- 
icant frequencies seen by Pinheiro et al. (2003) are not in 
disagreement with our results. 



4. The spectrum of V346 Ori 

4.1. Empirical determination of atmospheric parameters 

^From literature data we know that V346 Ori is classi- 
fied as a spectral type A5 or later and luminosity class 
III-V. For our classification we fo llowed the recipes de- 
scribed bv lGrav &: Garrisonl (|1989al lb[) for early F-type and 
late A-type stars, c o mplem ented with additional criteria by 
iJaschek fc Jaschekl (|l987f ). 

/From inspection of the spectral region typically used 
for classification (e.g. Fei at A4046 A; Srii at A4077 A; 
Feii-Tiii at AA4172-4179 A; Fei at A4187 and 4203 A; Fell 
at A4415 A/A4481 A (Mgii) we conclude that V346 Ori 
is of spectral type more similar to A9 than A7. Similarly, 
we tried to estimate the luminosity class of V346 Ori by 
analysing the morphology of the blends at AA4172-4179 
A and at AA4395-4400 A (Fe ii-Ti Il-Sc ll) with respect 
to A4203 A and A4271 A (Fei). We computed the ra- 
tio AA4172-4179/A4271 and AA4395-4400/A4271, which is 
about 2. This means that V346 O ri luminosity class is IV 
or V (see iGrav fc GarrisM3[l989al) . 

The following step was to determine more precisely the 
spectral type and simultaneously to estimate the rotational 
velocity (vsini) of V346 Ori. To this aim, we used the list 
of A-type stars with reliable and homogeneous vsini pub- 
lished bv [Rover et al.l (|2002|) coupled with high resolution 
spectra from the ELODIE archive (Moultaka al. 2004). In 
particular we compared (visually) the spectrum of V346 
Ori with several standard stars of spectral type A8 and 
A9 and luminosity class IV or V, as well as a variety of 
vsini. The best agreement was obtained with HD124675, 
an A8IV star with vsini~130 km/s and solar metallicity 
([Fossati et al.[ [2008[ : [Bush fc Hint3[2008t) . Figure [1 shows 
the superposition of V346 Ori and HD124675 spectra in 
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Fig. 5. Frequency analysis for V346 Ori calculated with periodOA adopting the weighting scheme of the data described in 
the text. The solid line shows the 99.9% significant level, corresponding to 12.6 times the noise level. In each panel, one 
peak (i.e. the labelled frequency) is selected and removed from the time series and a new spectrum is obtained. Then, 
all the derived periodicities are fitted simultaneously to refine frequencies, amplitudes and phases. The analysis of the 
residuals produces the subsequent periodogram. The bottom figure in the right panel displays the periodogram after the 
prewhitening with all the significant frequencies. 



the region used for spectral classification, as well as the dif- We estimate that V346 Ori is an A8IV(± 1 spectral 
ference between the two spectra. Figure [7] show the region class) star with vsini'^lSO km/s and approximately solar 
of Ha, which results partially filled and with an emission metallicity. 
in the core due to the presence of hot gas around the star, 

a typical feature of Herbig Ae stars. , This resuh agrees with some of the previous deter- 

minations, but disagrees with others. Considering spec- 
troscopic measurements, V346 Ori has been classified as 



S. Bernabei et al.: Multisite photometry of the pulsating Herbig Ae star V346 Ori 



1.5 



D 
X 

3 



1.0 



0.5 



' ' ' r 

MnlFel SrII H6 



T 1 1 r 



1 1 1 r 



Fell+Till Fel Cal 
Pel 




4000 



4100 



4200 



4300 



4400 



4500 



A(A) 
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A2IV (iMora et alJ 120011) . A5II Ie dHCTbiel llQgOD, A7III 
Gray k Cor ballv 199^), A7 III (iBlondel et al.l [200l . A8 



A 8V (jVieira et all l2003t) . and A9 

Hernandez etani2006r FI Thus, our estimate is consistent 



(iValenti et al.. .20001 



with the last three values, marginally with the A7III evalu- 
ation by Grav & Corballv (1998) and Blondel ct al. (2GG_6). 
It is inconsistent with iHerbid (|l960f ) and lMora et all ^OOU) . 
The origin of this disagreement in spectral types cannot be 
easily identified. We suggest two possible causes: 1) The 
modification of the spectrum due to non-photospheric ma- 
terial placed along the line of sight causing the star to 
appear earlier in some cases; and 2) The use of different 
lines/methods for the classification. For example, using Ca 
II K and/or H lines can often result i n earlier spectral 
type assignment for this type of stars (see lGrav fc GarrisonI 
11989a: .Grav fc Corballv 1998). Accordin g to the tempera- 
ture scale of [Kenvon fc HartmannI (|l995f ). our estimate for 
the spectral type of V346 Ori (A8±l spectral class) corre- 
sponds to Tefi=7580±250 K. Similarly, our estimated lumi- 
nosi ty class IV-V co rrespond approximately to log 3=3.5- 
4.3 (jSchmidt - Kaledil982i) . 



4.2. Spectral synthesis determination of atmospheric 
parameters 

To validate our spectral classification and to determine the 
physical parameters of V346 Ori, we have used synthetic 
spectra. 



Note that some authors do not report the "e" after the spec- 
tral classification as it is usually done. We verified that this is 
not due to a lack of detection of emis sion (in Ha) , b ut to vari- 
ous circumstances: e .g. the results by IValenti et al.l (|2000l ) and 
IBlondel et all (|2006l ) come from the UV etc. 



The approach we used in this paper was to minimize 
the difference between observed and synthetic H^, H-y and 
profiles. As goodness-of-fit test we used the parameter: 



1 V^/^obs 



(1) 



where N is the total number of points, Jobs and Ith are 
the intensities of the observed and computed profiles, re- 
spectively, and (5/obs is the photon noise. The synthetic 
spectra were generated in three steps: first, we computed 
the stellar atm osphere model by using the ATLAS9 code 
(jKurucd Il993f ) t hen, the stellar spectru m was synthesized 
using SYNTHE (jKurucz fc Avrettlll981[ ) and finally, the in- 
strumental and rotational convolutions were applied. The 
ATLAS9 code includes the metal opacity by means of dis- 
tribution functions (ODF) that are tabulated for multiples 
of the solar metallicity and for various microturbulence ve- 
locities. 

In this study we followed the method described in 
ICatanzarol ([2008') , since the intersection of the three 
iso-surfaces is expected to improve the final solution. To 
decrease the number of parameters, first we computed the 
Ve sin i of the star. To derive the rotational velocity, we used 
SYNTHE, using a guess model computed with the TcS and 
logg estimated in the previous sub-section, to reproduce 
the profile of Mgii A4481 A line and obtained a value of 
125 ± 10 km s^"'^, in agreement with our previous estimate 
of We sin i ^ 130 km s~^ . For this star the best fit has 
been obtained for the model computed with solar ODF for 
Teff = 7550 ± 200 K and log g = 3.5 ± 0.4. 

Hq, is also present in our spectral range, we did not 
use it principally for the evident emission it shows in the 
core. The synthetic line overplotted onto the observed 
one is shown in Fig. [51 It has been computed by using 
the atmospheric parameters derived from the other three 
Balmer lines. Also shown in the figure for comparison 
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Fig. 7. As in Fig. \6\ but in the region of Hq,. Note the 
off-center emission present in the difference spectrum. 

are the synthetic spectra with Toff ± 5 T and log g ± S log g. 

In conclusion, we have evaluated the atmospheric pa- 
rameters of V346 Ori by using three different methods: I) 
Stromgren-Crawford photometry; 2) empirical determina- 
tion of spectral and luminosity class; and 3) comparison 
with synthetic spectra. These approaches give consistent 
results within the errors. Given that the comparison with 
synthetic spectra is the most precise and physically based 
among the methods we used, in the following we will use the 
atmospheric parameters derived with this approach, but en- 
larging the error box in Tcff to ± 250 K to take into account 
the results of methods 1) and 2). 

5. Comparison with models 

Having evaluated the physical parameters of V346 Ori, 
we now proceed to constrain its mass and effective tem- 
perature (and, in turn, its luminosity) comparing the ob- 
served frequencies with model predictions. In order to re- 
duce the number of parameters in the theory versus ob- 
servations comparison, we decided to restrict the analysis 
to the empirically estimated ranges of luminosity (0.8 ^ 
logL/LQ <, 1.41) and effective temperature (7300-7800 K) 
for V346 Ori. These ranges have been derived on the basis 
of the atmospheric parameters estimated in the previous 
section and are represented in Fig. [9] by the dot-dashed 
box, together with th e evolutiona ry tracks computed us- 
ing the CESAM code (|More]||l997[) . From the intersection 
between the box and the evolutionary tracks we find that 
the predicted value of the stellar mass lies in the range 
from ~1.6 Mq to ~2.2 Mq. For each stellar mass in this 
range, we selected three stellar models corresponding to 
three different effective temperatures, separated by ~200 K, 
along the corresponding PMS evolutionary track, and per- 
formed the pulsation analysis in both radial and non-radial 
modes (1=0,1,2) using the ADIPLS cod^. Additional mod- 

http://astro.pliys.au.dk/ jcd/adipack.n/ 
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Fig. 8. Zoom of V346 Ori's spectrum in the region of H^, 
H^, H^ and Hq, in comparison with the best-fitting syn- 
thetic spectrum (Tcff = 7550 K, red line). Also shown in 
the figure for comparison are the synthetic spectra with 
Toff ± ST and logg ± S log 5 (blue dotted lines). Note 
that only H5, H^ and H^ were used for estimating the at- 
mospheric parameter of V346 Ori; Hq, is shown for compar- 
ison. 
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Fig. 9. HR diagram showing PMS evolutionary tracks com- 
puted with the CESAM code for the labelled stellar masses 
(in Mq). The dotted lines re present the large separation 
fit obtained bv iRuoppo et al.| |2007). The rectangular box 
delimits the range of luminosity and effective temperature 
of V346 Ori from the literature. The best fit model is given 
by the filled dot. 
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Fig. 10. Predicted frequencies as a function of the radial 
order n. Symbols are connected by a full line in the case 
of /=0, by a dashed line for 1=1 and by a dotted line for 
1=2. The horizontal lines represent the six most significant 
frequencies. 

els along the 2.3 M0 have been computed in order to ex- 
plore the behaviour of the predicted frequencies at higher 
stellar masses. The predicted frequencies, for each selected 
stellar model, are used to reproduce the six most signifi- 
cant periodicities (from /i to /g). The comparison between 
theory and observations is shown in Fig. 1101 where the fre- 
quencies are plotted as a function of the radial order for 
three assumptions 1=0,1,2, varying the stellar mass and 
the effective temperature. Inspection of this plot, together 
with the search for the minimum suggest that the 
best agreement is obtained for M=2.1 M©, logL/L© =1.37 
Teff =7314 K. The position of the best fit model in the HR 
diagram is represented by the large dot in Fig. [9l The de- 
rived value of the surface gravity of this model is log .g=3.8, 
close to our empirical estimate of log5=3.9±0.4. The vari- 
ous panels of Fig. [TU] show that increasing or decreasing the 
stellar mass and/or the effective temperatures cause the dis- 
crepancies between observed and predicted frequencies to 
increase. We have also verified that no reasonable match of 
the observed frequencies is found for stellar masses smaller 
than M=1.9 Mq. 

The comparison between the predicted and observed fre- 

* Given the observed frequencies of V346 Ori, we look for mod- 
els, into the empirical box of Fig. [O] having < 1- ^^s usual, 

the is defined as = -I7 ''/?>'°~'^t^j where A*' is the total 

number of observed frequencies, /obs and /th are the observed 
and the theoretical frequencies respectively, 5 fobs, Sfth are the 
corresponding error. If the program does not find any models 
with ^ li that means that not all frequencies are matched 
whithin the errors, it discards one by one the non-matched fre- 
quencies, computes again the for the remaining frequencies 
and searches again for models with x'^ ^ 1- The model with the 
least value of the is the best-fitting model. 



quencies in the echelle diagram is presented in Appendix 
B, confirming that we find a combination of mass, effective 
temperature and luminosity that is able to best reproduce 
the observed frequencies. We emphasize that this is true 
only within our explored range of stellar parameter and 
with the resolution in mass and effective temperature of 
the current model set 

6. Final remarks 

We have presented the results of a multisite campaign on 
the pulsating Herbig Ae star V346 Ori. The frequency 
analysis has allowed us to extract six highly significant 
frequencies on the basis of a no-weight Fourier analysis, 
and additional seven less significant frequencies by adopt- 
ing a "deviation weight" scheme. The measured period- 
icities have been compared with the predictions of non- 
radial adiabatic models. We found that the resulting best- 
fit model has a stellar mass of 2.1±0.2 M©, luminosity 
logL/L© = 1.37^q }3, log5=3.8±0.2, and effective temper- 
ature of ~7300±200 K, where the errors are upper limits 
based on the adopted mass and effective temperature steps. 
These estimates are in agreement with the spectroscopic es- 
timates, within the associated uncertainties. 

The value of the asteroseismic luminosity is signifi- 
cantly higher than that obtained by other authors. For 
example, Hernandez et al. (2005, H05 hereafter) derived 
a value of logL/L© = 0.88 (no errorbar is given) using 
the average astrometric Hipparcos distance (~330± 15 pc) 
to the O ri OBla association a nd the PMS evolutionary 
tracks bv lPalla fc Stahleil (|1993D . assuming Av=0 for V346 
Ori. However, the uncertainties involved in this calculation 
are large. Indeed, to derive the luminosity of a star, even 
knowing the distance, one should know accurately the vi- 
sual magnitude and the interstellar extinction. This is not 
straightforward for V346 Ori, which sh ows a wide var i- 
ability in the visual band (V=9. 8-11.8: iThe et al]ll994l ). 
due to the presence of circumstellar material around the 
star (see e.g. the SED published bv Hernandez et al. 20061). 
Simil arly the value of 4^ ^® ^sed is not easily achiev- 
able: iHernandez et al.l (|2005[ ) used Av=0, but the aver- 
age absorption for the sub-group 25 Ori to wh i ch V 346 
Ori is beheved to belong to (see iBricefio et al.l l2007l ) is 
< Av >= 0.29. It is useful to recalculate the logL/L© 
value for V346 Ori adopting the distance by H05, but tak- 
ing into account the value at minimum light for V346 Ori 
V=9.8 (which should be the closest to the photosphcric 
one), and a visual abs orptio n Ava=0.41±0.1 as provided by 
the maps of Schlegel et al.l (|l998l ). which arc valid in an 
area of 6.1'x 6.1' around the star. As a result, allowing for 
an error of 5% for the distance and 10% for both the visual 
magnitude and absorption, we get logL/L0=1.17±O.ll, 
which agrees to within Icr with the results of our astero- 
seismological analysis. Of course, a better agreement could 
be obtained either by placing V346 Ori slightly beyond 
the Orion OBla association (indeed the distance implied 
by the present results is D=400±50 pc that corresponds 
to the revised distance of the Orion Nebula Cluster, e.g. 
Jeffries, 2007) or supposing that the extinction due to cir- 
cumstellar material has been underestimated. In the first 
case V346 Ori would not be a member of the Orion OBla 
association. To our knowledge, the membership of V346 Ori 
to OBla has not been confirmed spectroscopically. To cast 
light on this complicated scenario, an accurate measure- 
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ment of the radial velocity of V346 Ori would be highly 

desirable. 
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Table A.l. Log of the observations. 



HJD-2450000 
start (days) 


HJD-2453000 
end (days) 


Duration 
(hours) 


Telescope 


677.657 


677.752 


2.3 


REM 


678.654 


678.721 


1.6 


REM 


683.641 


683.734 


2.2 


REM 


684.638 


684.769 


3.2 


REM 


685.635 


685.779 


3.4 


REM 


689.690 


689.751 


1.4 


REM 


690.622 


690.734 


2.7 


REM 


695.421 


695.486 


1.6 


Loiano 


695.608 


695.791 


4.4 


REM 


696.605 


696.730 


3.0 


REM 


698.626 


698.744 


2.8 


REM 


699.626 


699.741 


2.8 


REM 


700.626 


700.738 


2.7 


REM 


701.598 


701.723 


3.0 


REM 


702.604 


702.733 


3.1 


REM 


703.395 


703.472 


1.8 


Loiano 


703.589 


703.769 


4.3 


REM 


704.587 


704.675 


2.1 


REM 


704.724 


705.032 


7.4 


SARA 


705.597 


705.750 


3.7 


REM 


707.578 


707.710 


3.1 


REM 


708.605 


708.717 


2.7 


REM 


709.572 


709.714 


3.4 


REM 


710.585 


710.711 
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Appendix A: Log of the observations 

Appendix B: Theory versus observations in the 
echelle diagram 

A powerful tool to investigate the prediction capability of 
pulsation models is represented by the echelle diagram, 
where frequencies are reduced modulo A/ by expressing 
them as fn,i = fo + A/ + / where A/ is the large separa- 
tion, fo is a suitably chosen reference, and k is an integer 
such that / is between and A/. The echelle diagram in- 
cludes the information on the frequency large separation so 
that the comparison between modeled and observed peri- 
odicities in this diagram can allow us to constrain the prop- 
erties of the frequency spectrum. Hence, for all the mod- 
els located in the empirical box of Fig[51 we have used the 
computed frequencies to build a theoretical echelle diagram 
and search for the one that best reproduces the observed 
frequencies. As a result we obtain the best fit model quoted 
in Sect. 9. Looking at the behaviour of the best fit model 
in Fig. [To] we notice that the most significant frequency /i 
can be reproduced within ± 5/iHz cither with 1=0 or with 



1=2. Assuming 1=0, as also supported by the idea that the 
largest amplitude frequency (/i) is a radial mode, we obtain 
the echelle diagram shown in the left upper panel of Fig. IB. II 
for a large frequency separation A/ =35.8 /iHz. This quan- 
tity is evaluated using the analytical relation between the 
large frequency separation, the s tellar luminosi t y and the 
effective temperature derived bv iRuoppo et al.l (|2007D . As 
already found in Fig.fTUl the best fit model is able to repro- 
duce the most significant frequencies within an uncertainty 
smaller than 5 fj, Hz apart from f^. This frequency is ver y 
close to /s and, as discussed bv iBreger &: Bischol (|2002[ ). 
close frequency pairs should be investigated accurately in 
order to establish if they are real or not. 

As already specified, the best fit pulsation model has 
been found by requiring it to match only the six most sig- 
nificant frequencies. If the other seven measured frequencies 
(from /7 to /13) are added in the echelle diagram we ob- 
tain the plot shown in the left lower panel, where the only 
deviating frequency, beyond the already noted /4, is /12. 

From our spectroscopic analysis we have found that 
V346 Ori has a projected rotational velocity of v sini = 
125±25 km/s. The effect of rotation on non-radial modes 
can be roughly estimated using the asymptotic relation for 
non-radial displacements. This is given by mfrot, where m 
is the spherical harmonic number and frot = 2TiR's^ini 
rotational frequency. For the inclination angle i we consid- 
ered the two cases ii = ^ and 12 — ^, while for the radius 
we took the value of our best fit model. Then, with these as- 
sumptions, the resulting rotational splitting is A/i ~ m5.3 
fiHz and A/2 ~ m6.2 /iHz, for vsini =100 km/s and 
A/i ~ m7.9 /iHz and A/2 ~ rn9.3 /xHz for vsini =150 
km/s. It is interesting to note that the displacement of /4 
and /12, from the I = 1 sequence seems to be consitent 
with the rotational splitting estimate for vsini =150 km/s, 
with these frequencies corresponding to Z = 1, m = ±1 
modes. However, we are aware that this is a very simpified 
treatment of the rotational splitting phenomenon and that 
a more accurate analysis also by means of asteroseismologi- 
cal models including rotation is needed in order to properly 
test this possibility. 

If we discard the assumption that /i corresponds to a ra- 
dial mode and explore the alternative possibility suggested 
by the behaviour of the best fit model in Fig. (TUl i.e. that 
/i is a 1=2 mode, we obtain the plot shown in the right 
upper panel of Fig. IB. II In this case all six frequencies are 
reproduced within ±5 /iHz, apart from /a, that deviates a 
bit more from the 1=2 sequence. However, if the other seven 
less significant frequencies are added, we find that theoret- 
ical predictions are not able to provide a good match for 
all the frequencies (see the right lower panel of the same 
figure) . 
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Fig. B.l. Echelle diagram for the model M=2.1 M©, Teff=7314 K compared with the frequency data. The parameters 
adopted to match the data are labelled. Left: with the assumption that the largest amplitude frequency (/i) is radial 
(1=0) and by matching only the six most significant frequencies (upper panel) or all the frequencies (bottom panel). 
Right: with the assumption that the largest amplitude frequency (/i) is non radial (1=2) and by matching only the six 
most significant frequencies (upper panel) or all the frequencies (bottom panel). 



